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Abstract 

Using a mixture of Unear algebra and statics, we derive what can be viewed as a 
slight generalization of the Coppersmith-Tetali-Winkler Identity H(i,j) + H(j,k) + 
H(k, i) = H{j, i) + H(k,j) + H{i, k) for hitting times of a random walk. 

1 Linear Algebra 

In this paper, the bold-faced characters are reserved for denoting vectors; the dimension of 
the vectors is n > 3, except when the context suggests otherwise. For example, f is a vector 
with coordinates /i, . . . , /„, while and 1 are the all- zero and all-one vectors, respectively. 
Consider the system of linear equations 

(1) Ax = f 

where x is the vector of the unknowns, f is a fixed vector, and A is an n x n matrix such 
that: 

(i) A is symmetric; 

(ii) the off-diagonal entries of A are non-positive; 

(iii) A is irreducible, which is to say that the simple graph whose adjacency matrix has the 
same off-diagonal zero pattern as A is connected (throughout the paper, we refer to 
this graph as "t/ie underlying graph G"); 

(iv) Al = 0. 

Applied to a matrix with properties (i)-(iii), the Perron- Frobenius Theorem states that its 
least eigenvalue has multiplicity one, and the corresponding eigenvector can be chosen to 
have all its coordinates positive. Hence, by (iv), A is positive semidefinite and of rank n — 1. 
Consequently, (1) has no solution if f ■ 1 7^ and a solution unique up to a shift by a multiple 
of 1 if f ■ 1 = 0. For 2 e {1, . . . ,n}, set f ^ := (A, . . . , /, - f ■ 1, . . . , /„) (so that f . ■ 1 = 0) 
and let X j = {xn, . . . , Xj„) denote the [unique] solution to the system of linear equations 

(2) Ax=f, 
such that Xii = 0. Our goal is to demonstrate the identity 

(3) a;i2 + X23 + a;3i = a;2i + X32 + X13. 
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Suppose first that n > 3 and let us apply [symmetric] Gaussian elimination to (2), to 
eliminate the non-zero off-diagonal entries of A in the last column and row. The reader can 
easily see that the first n — 1 unknowns are not affected by such an elimination, while the 
upper left (n — 1) x (n — 1) corner, A', of the matrix obtained from A, and the (n — 1)- 
dimensional vector, i[, whose entries are the same as the first n — 1 entries of the vector 
obtained from f j, have the same properties as A and f j, respectively. In other words, the 
system of n — 1 linear equations 

(4) A'x'^f: 

with n — 1 unknowns be considered in lieu of (2), and if Xj = {xn, . . . , Xin) is 

a solution to (2) then := {xn, . . . , Xi^n-i) is a solution to (4). Hence, in order to establish 
(3), we may assume that n = 3. 

Still, to verify (3) for n = 3 directly is a rather tedious task. After all, one would have 
to solve three systems of three equations each in general. Instead, we prefer the following 
detour. In the next section, we will see that (3) is trivially true if the underlying graph G 
is the star Ki^n~i- If n = 3 then G is either the star Ki^2 or the triangle K^. In the latter 
case, the matrix A has the form 

Oi2 + Oi3 —as —a2 
—as ai + as — ai 
—a2 —ai «! -|- a2 

where ai, a2, and are some positive numbers. Set c := aia2 + a2as + asai and Pi :— c/ai, 
i e {1, 2, 3}. Then, the matrix 



B : = 
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P1 + P2 + P3 



has properties (i)-(iv) and it is not difficult to verify that A is the upper left 3x3 corner of 
the matrix obtained from B by symmetric Gaussian elimination. Hence, (3) holds for A if 
it holds for B. Which it does, since the underlying graph for B is the star Ki ^- 

2 Mechanical Model 

The systems (1) and (2) have the following mechanical interpretation: n masses labelled 1 
through n are submerged in water, so that each mass j is subject to a (positive or negative 
or zero) weight, fj. In addition, every two masses j and k are connected by a rubber band 
of resistance ajk = a^j > where —ajk is the jk entry of A. (A rubber band is of resistance 
at most a > if and only if a force of magnitude al is sufficient to stretch it to length I.) 
If mass i is nailed to the origin, such a system of masses has a unique equilibrium. Let Xij 
denote the altitude of mass j in this equilibrium. (So that Xu — 0.) The fact that the system 
is in an equilibrium is expressed by equating the sum of the forces acting on every mass j to 
zero. In an equilibrium, these forces are the weight /j, the reaction — f • 1 = — (/i -|- . . . -|- /„) 
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of the nail if j = i, and the resistance forces aij^xa — Xij), . . . , anj{xin — Xij). In other words, 
the vector (xji, . . . , is the unique solution x j to the system (2) such that Xa = 0. 

Using this physical interpretation of (2), we see that (3) is trivially true if the underlying 
graph G is the star Ki^n-i with say center k. This is because Xij = Xik + Xkj in this case. 
This completes the proof of (3). 

Remark 1: In the mechanical model, we could replace the weights f] by arbitrary d- 
dimensional forces and, correspondingly, the altitudes Xij by the ci-dimensional position 
vectors. Then, the "(i-dimensional generalization" of (3) would be proved coordinate-wise. 

Remark 2: The identity Xj^jj + Xjjig + . . . + Xi^i-^ = Xi^i^ + Xjgjj + . . . + Xi^i^, can be shown 
for any sequence ii, . . . ,ik of indices using (3); we leave this to the reader. 



3 Random Walks 

Let G be a connected graph on n > 3 vertices. To obtain a random walk out of vertex j in 
G, one recursively builds a sequence zoii«2 ... of vertices, in which io = j and, for t >0, the 
vertex i^+i is chosen among the neighbors of it with the uniform probability. The hitting 
time H{j, i) is the expectation of the least number t such that it = i. (So that H{i, i) = 0.) 
If 2 7^ j then, clearly, 

(5) HU,^) = l + ^ E H{k,^). 

Multiplying both sides of (5) by deg(j), we see that the vector x, := . . . ,H{n,i)) 

is a solution to the system of linear equations 

(6) Ax=g, 

where gj is some vector coinciding with the vector f := (deg(l) . . . , deg(ri)) in all the 
coordinates except perhaps i, and A is the difference between the diagonal matrix diagjf } 
and the adjacency matrix of G. Clearly, A satisfies (i)-(iv), whence gj must be the vector 
fj obtained from f as in Section 1. Hence (3), taking in this context the form of the 
Coppersmith- Tetali- Winkler Identity H{i,j) + H{j,k) + H{k,i) = H{j,i) + H{k,j) + 
H{i, k), holds. Also, substituting x , into equation i of (6) gives the identity 

J2 H{k,i) = 2m-degii) 

{k:ikeG} 

where 2m := deg(l) -|- . . . + deg{n) is twice the number of edges in G. It follows that the 
return time R{i), i.e. the expectation of the smallest natural t such that it = i in a random 
walk out of i, is 

1 /.NX 2m 

deg{i) deg{i) 
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